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Abstract 

This is our response to a comment by Walter Eifler on our paper 'A simple model 
for the short-time evolution of near-surface current and temperature profiles' 
(arXiv:physics/0503186, accepted for publication in Deep-Sea Research IT). Al- 
though Eifler raises genuine issues regarding our model's validity and applicability, 
we are nevertheless of the opinion that it is of value for the short-term evolution 
of the upper-ocean profiles of current and temperature. The fact that the effective 
eddy viscosity tends to infinity for infinite time under a steady wind stress may not 
be surprising. It can be interpreted as a vertical shift of the eddy viscosity profile 
and an increase in the size of the dominant turbulent eddies under the assumed 
conditions of small stratification and infinite water depth. 
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Eifler (2005) criticises our model ( Jenkins and Ward 200,4 referred to as JW) 



for a number of reasons: 

(1) If our model is run with a step- function wind and wind stress applied at 
time t = 0, the velocity field u(z,t) becomes 



u(z, t) = XU exp — - 
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and satisfies 



du u* 2 z d 2 u 

dt = ~Wd^' ( ' 

where z is the downward-pointing vertical coordinate, \U is the interface 
velocity, U is the wind speed, and u* is the waterside friction velocity. In 
JW we state that the quantity 

e -— (3) 
pp XU 

plays a role similar to that of an eddy viscosity. Eifler disagrees with 
our interpretation, using a stricter definition of eddy viscosity e, in our 
one-dimensional case it is defined using the equation 

du d ( du\ , . 

from which he deduces from (2) that 
u* 2 fu* 2 V 

so that the vertical momentum flux 'has to be sustained by an eddy 
viscosity tending towards an infinite value'. 

(2) Eifler computes the time dependence of the mixed-layer depth (defined 
for the JW model as the depth where the local vertical momentum flux is 
a specified fraction of the surface momentum flux), and compares it with 
an 'empirically-tested standard formulation'. The mixed-layer depth z m i 
of the JW model obeys 

j^T=x% 

Zml = 16 U*t, (6) 

where A T=X% is defined by Eifler's Eq. 10. The empirically-tested stan- 
dard formulation gives a mixed-layer depth D m \ which obeys 

D ml = 10.5 M *t 1/2 , (7) 

the units of D mi , u*, and t being metres, ms _1 , and seconds, respectively. 
The different time dependence of the mixed layer depths in Eqs. 6 and 7 
is stated to be a reason against the validity of the JW model. 

(3) The derived numerical thermal model of JW shows in Fig. 3 (of JW) 
a mixed layer depth which develops linearly in time, and the agreement 
with observations is not particularly good in this case, giving further 
evidence against the usefulness of the JW model. 
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Although Eifler raises genuine issues which question the validity and appli- 
cability of the JW model, we are nevertheless of the opinion that it has a 
degree of usefulness: indeed, Eifler does state in his conclusion that JW 'have 
provided a promising analytical solution of the vertical transport equation . . . 
with very interesting mathematical properties'. 

Regarding Eifler's Reason 3 above, we concede that the model (not surpris- 
ingly, given its simplicity) does not always give good agreement with the mea- 
sured data, although the particular reasons for this may also include horizontal 
advection effects. Of course, perfect dynamical agreement cannot in any case 
be obtained, since the JW model assumes zero Coriolis force, neutral stratifi- 
cation, and infinite water depth. 

Reason 1, that the eddy viscosity e tends to infinity for infinite time under a 
steady applied wind stress, is not such a serious objection as it first appears. 
Equation 5 may be re-written as 

u" 2 r fu' 2 \ 



rom which we see that the graph of e versus z shifts upwards at a velocity 
u* 2 /(XU) . This result may also be interpreted as an increase in the length 
scale of the dominant turbulent eddies which impact the surface (and other 
depths). The increase of e to infinity at infinite time is a consequence of the 
depth being infinite. The fact that the near-surface values of the eddy viscosity 
are no longer close to linearly proportional to depth is in fact reminiscent of 
the breaking- wave enhance d values o f near- surface eddy viscosity shown in 
Figs. 3.18, 3.19, and 7.16 of lBurchardl (|2002h . 



Reason 2, that the JW linear t dependence of mixed layer depth z m i is in dis- 
agreement with the empirically-tested t 1 / 2 dependence of mixed la yer depth 
Dmi, is rather interesting. Eq uation 7 cor responds to Eq. 6.1 in iBurchard 



(2002), which was sugg estedbv lPricel(|l979h . The constant 10.5 in (7) is in fact 



I.OBNq 1 ^ 2 , with N Q , the initial Brunt- Vaisala frequency, being set to 10~ 2 s _1 . 
For an initially linearly-stratified water column, or, in fact, for a water column 
with any stably-stratified density profile, subjected to a constant input of tur- 
bulent energy at the surface for mixing, the result D m i cx t 1 / 2 is a consequence 
of the requirement that the potential energy of the water column increases at 
a constant rate due to mixing. Similarly, the result z mi oc t is a consequence of 
a constant rate of momentum input to the water column. If the Coriolis force 
due to the Earth's rotation is taken into account, z m i can be expected not 
to increase indefinitely, as the wind-induced current would be restricted to a 
boundary laye r of Ekman depth typically of the order of magnitude of u*/f 



f)Madsenl 119771 ). Further study is required to resolve these different boundary 



layer time dependences, and to evaluate the joint effects of stratification and 
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rotation within the simple JW model framework. 
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